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In this paper we apply the self-consistent generalized Langevin equation theory (SCGLE) of 
dynamic arrest for colloidal mixtures to predict the glass transition of a colloidal fluid permeating 
a porous matrix of obstacles with random distribution. We obtained the transition diagrams for 
different size asymmetries and so we give an asserted description of recent simulations results [K. 
Kim, K. Miyazaki, and S. Saito, Europhys. Lett. 88, 36002 (2009)] of Quenched- Annealed and 
Equilibrated-Mixture systems which reveal very different qualitative scenarios which are in apparent 
contradiction with theoretical predictions of Mode Coupling Theory (MCT) [V. Krakoviack. Phys. 
Rev. E 75, 031503 (2007)]. We show that SCGLE theory predicts the existence of a reentrant 
region in EM systems as predicted using MC theory. However, opposite to MCT predictions, we 
show that it is practically impossible to distinguish a rentrant region in QA systems if it would 
exist. Qualitative comparisons are in good agreement with simulation results and thus, we propose 
SCGLE theory as a useful tool for the interpretation of the arrest transition in ideal porous systems. 

PACS numbers: 23.23.+X, 56.65.Dy 
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The study of static and dynamic properties of fluids 
permeating porous materials is an interesting topic that 
has attracted increasing attention due to both its rel- 
evance in the study of systems in biology, chemistry, 
physics and engineering [l[ and its scientific importance 
in the understanding of diffusing phenomena [2 H| . For 
this, the development of theoretical schemes describing 
this class of systems has deserved considerable work from 
many perspectives. In this sense, a binary mixture in 
which one of the species remains immobile while the sec- 
ond species diffuse constitutes a model system that has 
been widely employed not only in theoretical or simula- 



nomena as well as for its potential applications in in- 
dustry. Nevertheless it has attracted a lot of interest 
in the last decades 12h14|. there are still many open 
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tion studies [J-[9| but also in experimental works [1 
as an idealization of a fluid immerse in a porous medium. 
Despite the simplicity of this model, some parameters like 
size asymmetry or matrix structure can be selected in or- 
der to mimic different systems. There are two protocols 
widely used to generate the obstacle matrix in order to 
simulate different classes of disordered porous materials, 
known as Quenched- Annealed (QA) Mixture and Equi- 
librated Mixture (EM) [8]. In the first case, the porous 
matrix is prepared by quenching spherical particles of a 
monocomponent system in a typical equilibrium config- 
uration and after that, mobile particles are inserted into 
the void spaces of the matrix of immobilized particles. 
By the opposite, in the EM protocol, all particles equili- 
brate together and then, a fraction of them arc suddenly 
immobilized and so, the structure of the obstacle matrix 
corresponds to that of one species of a mixture in equi- 
librium. 



On the other side, the study of the structural glass 
transition of fluids is another topic of interest because 
of its importance in the comprehension of dynamic phe- 



problems regarding to the nature of this dynamic phe- 
nomenon. For this, it is always desirable to have theoret- 
ical schemes that allow for the interpretation of experi- 
mental and computer simulation results. In this sense, 
Mode Coupling Theory (MCT) [ill - fl7t has shown its ca- 
pability in predicting liquid-glass transitions in a vari- 
ety of systems and conditions [l8l - |20| . In this context, 
theoretical studies of liquid-glass transition of fluids con- 
fined in porous media deserves special interest in order 
to achieve a better understanding of the slow dynamics 
of these systems. Recently, based on the replica Orstein- 
Zernike equations [i , 21 , 22 1 , Vincent Krakoviack has ex- 
tended the Mode Coupling Theory to consider the glass 
transition of a fluid immersed in a porous medium Q. 
Within this Replica Mode Coupling Theory (RMCT) he 
obtained the dynamic arrest diagrams for QA mixtures 
and showed that RMCT predicts a reentry phenomenon 
for matrix densities higher than the localization thresh- 
old. This contribution represents a crucial first step in 
the theoretical description of the slow dynamics of a fluid 
immersed in a porous medium. 

Because of the complexity of the problem and the va- 
riety of scenarios regarding the glass transition in porous 
media, computer simulation studies have served to re- 
veal some special features that allow to evaluate theo- 
retical predictions. In particular, recent work of K. Kim 
and coworkers Q has shown the strong dependence of 
the glass transition diagram on the protocol in which the 
porous matrix is prepared. In their work, they performed 
simulations of hard sphere systems permeating an obsta- 
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cle matrix prepared using both QA and EM protocols, 
and compared their results with RMCT predictions. Two 
important aspects should be addressed: a) by construc- 
tion, there is no way of predicting EM glass transition 
diagrams using RMCT and b) simulation results are in 
apparent contradiction with RMCT prediction about the 
reentrant behaviour in QA systems. 

In this sense, the prediction of dynamic arrest tran- 
sition diagrams employing an alternative theoretical 
scheme is crucial for the interpretation of the simulation 
results mentioned above. The Self Consistent General- 
ized Langevin Equation theory (SCGLE) has been de- 
veloped in the context of colloidal systems fisT - E^ and 
it has shown their capability in describing the dynamic 
properties of a variety of systems [26j]. The extension of 
this theory to mixtures has been carried out [2^, HI] and 
allows us to consider QA an EM systems, just assuming 
that one species does not diffuse. In fact, it has already 
been employed to describe the time evolution of the col- 
lective diffusion properties of a fluid immersed in an ideal 
porous medium using both protocols [f|. Moreover, there 
have been derived numerical criteria for the prediction of 
the liqu id-glas s transition for several systems and con- 
ditions 29|432|, requiring as input only structural prop- 
erties as the static structure factors, S a p(k). Actually, 
it has been shown that, for some cases, this theory is as 
useful as MCT in the localization of dynamic arrest tran- 
sitions. Since SCGLE theory has been applied succesfully 
in the description of the dynamics of colloidal fluids per- 
meating ideal porous media as well as in the prediction of 
the liquid-glass transition in mixtures, it is obvious that 
the next step in this direction is the extension of this the- 
ory to predict the dynamic arrest transitions of QA and 
EM systems. This is the aim of this work. 

The relevant dynamic information of an equilibrium v- 
component colloidal suspension is contained in the v x v 
matrix F(k,t) whose elements are the partial intermedi- 
ate scattering functions F a p(k,t) = (n a (h, t)np{— k', 0)) 
where n a (k, t) = ex P[*k-r i (i)]/y / 7Va, with r^t) be- 

ing the position of particle i of species a at time t. The 
initial value F a /s(k, 0) is the partial static structure factor 

s a p(k) EE 111 . 

The multi-component version of the SCGLE the- 
ory consists of a set of exact time-evolution equations 
that governs the relaxation of the partial intermediate 
scattering functions F a p(k,t) and their self contrapart 
F$(k,t) = 6 a p (exp [tk • AR( Q )(f)]), where AK^(t) is 
the displacement of any of the N a particles of species a 
over a time t, and 5 a p is Kronecker's delta function. The 
Laplace transform (LT) F(k, z) of the matrix F(k, t), can 
be written as \2l\ , 

F(k, z) = {z + (I + C(k, z)Y 1 k 2 DS- i y 1 S, (1) 

where the elements of the matrix D are given by D a p = 
SapD^, with D° being the diffusion coefficient of species 



a in the absence of interactions. This is related with 
the solvent friction coefficient on an isolated particle 
of species a, (®, through the Einstein relation, D° = 
ksT/C,^. The elements C a p{k,z) of the matrix C{z) 
are the LT of the so-called irreducible memory func- 
tions C a p{k,t) 30]. The corresponding result for the 

"self" component, F^(k,t). Defined as F^(k,t) = 
<W (exp [ik • AR< Q '(i)]), where AR( Q )(t) is the displace- 
ment of any of the N a particles of species a over a time 
t, and S a /3 is Kronecker's delta, in Laplace space is 



F^{k,z) = {z+(I + C^(k,z))- 1 k 2 D} \ (2) 

where the matrix C^ s '(Jk, z) is the corresponding irre- 
ducible memory function which is related with C(k,z) 
by 



C(k,z) = C* (s) (fc,z) = A(fc)AC*(z), 



(3) 



The matrix A£*(z) in is a diagonal matrix whose a- 
th diagonal element, A£*(z), is the Laplace Transform 
of the time-dependent friction function of particles of 
species a. Such function reads [3l| 



D° 
3(2tt) s 



d 3 kk 2 [F^ (k, t)] aa [cy/n~F(k, t)S^V^h] c 

(4) 



with the elements of the fc-dependent matrices h and 
c being the Fourier transforms h a p(k) and c a p(k) of the 
Ornstein-Zernike total and direct correlation functions, 
respectively. Thus, h and c are related to S by S — 
I+^/nh^/n = [I — ^/nc^/n] -1 , with the matrix y/n defined 
as [\/n\ a fi = 5 af3 ^/n^. 

As illustrated in Ref. SH] , the solution of the SCGLE 
theory provides the time and wave- vector dependence of 
the dynamic properties of the system contained in F(k, t) 
and F^(k,t). It also provides equations for their long- 
time asymptotic values, referred to as non-ergodicity pa- 
rameters, which play the role of order parameters for the 
ergodic-non-ergodic transitions. The most fundamental 
of these results 3l) is the equation for the asymptotic 
mean squared displacement 7 Q = lim^oo < JAR' ') 2 >, 
which reads, 



3{2n)3J d°kk»{X[A + k^} aa 



(5) 



where S is the matrix of partial static structure factors, 
h and c are the Ornstein-Zernike matrices of total and 
direct correlation functions, respectively, related to S by 
S = I + y/nhy/n = [I — *Jncy/n\~ l , with the matrix 
yjn defined as [v^c^ = 5 a p-\/na, and A(fc) is a diagonal 
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matrix given by \ a p{k) = 5 a p[l + (fc/fce ) 2 ] 1 , where 
is the location of the first minimum following the 
main peak of S aa (k). 

To generalize the SCGLE theory to calculate the dy- 
namic properties F a p(k, t) of a colloidal fluid mixture of 
/j, components immerzed in a porous matrix, we use the 
same model employed in Ref. @. This consists of a 
j/-component colloidal mixture where a fraction of the 
particles, /i of the v species, are diffusing in the random 
matrix of obstacles formed by restant v — /a immobile 
(self-diffusion coefficients identically zero) especies which 
play the roll of porous medium. The derivation of corre- 
sponding equations to calculate the partial intermediate 
scattering functions F(k,t), F < -s)(k,t) and the long time 
mean squared displacement j a = limt-^oo < (AR^) 2 > 
of the mobile species is inmediate. Simply we make the 
self-diffusion coefficients = of all species that rep- 
resents the immobile particles in the system of Eqs. 11121 
and |H the result is a reduced set of \i equations for the 
fj, mobile species, with the same functional form. 

We will consider a binary hard sphere mixture of N\ 
mobile particles of diameter a\ and N2 obstacle parti- 
cles of diameter 02. The only macroscopic control pa- 
rameters are the volume fractions <j>i = irriicrf/6, where 
rii = Ni/N and N — N\ + N%, and the size asimme- 
try defined as 5 — a\l<Ji. The information about the 
mobility of the particles is contained in the self diffusion 
coefficients, so that, making D\ = in Eq. djimplies 
that particles of species 2 remain immobile. On the other 
hand, the structure of the porous medium is contained in 
the matrix of static structure factors, S. In practice, the 
selection of the protocol to generate the porous matrix 
determines the way on calculating the elements of matrix 
S. In this work, we used replica Orstein-Zernike equa- 
tions with given as the input to obtain the rest of the 
elements of the matrix S. Thus, for SmQs) we employed 
the static structure factor of a monodisperse system in 
the case of QA mixtures, and a static structure factor of 
a binary mixture in the case of EM systems. 

In Fig. [I] we compare the dynamic arrest diagrams ob- 
tained using Eq. |6]with both protocols to explore the sen- 
sitivity of our theory to the morphology of the medium. 
Dashed lines correspond to Quenched-Annealed (QA) 
and solid lines correspond to Equilibrated-Mixture (EM) . 
Three size asymmetries were studied: 5 = 0.5, 1 and 2; 
upper (black), middle (red) and lower (green) curves re- 
spectively. Observe that EM protocol shows an evident 
reentrance while for QA mixtures it is practically impos- 
sible to distinguish the existence of a reentrant region. 
We argue that this difference arises from the fact that 
thermalization process in EM systems, in which both 
species equilibrate together, leads to a self-generation of 
free- volume by mobile particles. This is easy to under- 
stand considering the fact that an increase in the number 
of fluid particles for a fixed number of obstacles implies 



an increase in the total density which means less volume 
occupied by particle and so, in certain region of the tran- 
sition diagram, once the obstacle particles are quenched, 
the fluid has more free volume available to diffuse. On 
the other hand, in QA systems, the volume occupied by 
the obstacles is the same independent of the concentra- 
tion of fluid particles. So, an increase on the number of 
fluid particles would not increase their mobility, in prin- 
ciple, and for that reason there is not a reentrant region 
for QA systems. These predictions are in agreement with 
comput" ^ : - @ 
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FIG. 1: Dynamic arrest phase diagrams of a monodisperse 
Hard Sphere fluid immersed in a random porous medium pre- 
pared in two diferent models: Quenched-Annealed (dashed 
lines) and Equilibrated-Mixture (continuous lines). And three 
size asymetries ratios between fluid and matrix particles: bot- 
tom S = 0.5 (black), intermediate 5 = 1 (red) and top 5 = 2 
(green). (Color on line) 



Qualitative comparisons between SCGLE predictions 
and computer simulation results are shown in figure [5] 
where we plot the SCGLE transition line (continuous 
line) and the fluid (solid) and arrested (open) states de- 
termined by computer simulation experiments of Kim 
et. al. for QA mixtures (upper panel) and EM systems 
(lower panel). Vertical axis corresponds to the obstacle 
volume fraction scaled with the percolation volume frac- 
tion, <pp, and horizontal axis corresponds to the mobile 
particles volume fraction scaled with the monocompo- 
nent glass transition volume fraction, <p p . For reference, 
we also plotted (dash line) the arrest transition line pre- 
dicted by Krakoviack with RMCT for QA systems. It de- 
viates considerably from simulation results and predicts 
a reentrancy that cannot be appreciated in simulation ex- 
periments. He explained this reentrancy in terms of the 
derealization of mobile particles caused by occasional 
collisions with other fluid particles trapped in neighboor 
cages. We believe that this effect -if there exist- could 
not be responsible of the reentrant pocket, as it has been 
pointed out in Ref. [8| where the authors have shown that 
it is hard to identify the existence of a reentrant region 
in the arrest transition diagram obtained with molecular 
simulations. 
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FIG. 2: Qualitative comparison between dynamic arrest di- 
agram of hard sphere fluid absorbed in a random porous 
medium calculated with SCGLE theory (continuous line) and 
simulations of Kim et al (fill symbols fluid phase and empty 
symbols arrested phase), for both QA (upper panel) and EM 
(lower panel) systems. RMCT prediction (dash line) is shown 
by reference. 



To understand the reentrant behavior of arrest transi- 
tion diagrams with both protocols, one must carefully an- 
alyze the mechanisms that give rise to this phenomenon 
in both EM and QA systems. First, Krakoviack proposed 
a mechanism to explain the notorious reentrancy in the 
diagram predicted by RMCT with QA protocol [7J , argu- 
ing that the dephasing of trajectories of the caged mobile 
particles by collision with other mobile particles causes 
an increase in the percolation threshold. Although this 
mechanism could explain the existence of a reentrancy, 
says nothing about the magnitude of this effect. We 
think, in agreement with arguments of Kim et. al. Q, 
that even when those collisions can eventually occur, the 
structure of the porous matrix is not modified and con- 
sequently the effects on the phase diagram should not 
be very noticeable. Thus, if this reentrant region does 
exist in the arrest diagram, the mesh size of the simula- 
tion should be finer to show it clearly. This fact could be 
observed in Fig. [5] where the transition line predicted us- 
ing SCGLE theory shows a barely noticeable reentrance. 
On the other hand, in the case of the EM protocol, the 
reentrancy is well_ understood and has already been ex- 
plained in Ref. Q. The fact that both species (mobile 
and immobile) are equilibrated together before stopping 
the obstacle particles, implies that increasing the concen- 



tration of mobile particles, due to purely entropic effects, 
they generate their own space to continue moving. For 
that, unlike the case of QA, in EM systems the structure 
of the porous matrix is affected by this increase and so, 
the reentrance is hightly visible, as shown by MD results 
and theoretical (SCGLE) curve. 

In summary, we have applied SCGLE to predict the 
dynamic arrest transition in a fluid permeating a porous 
matrix using Quenched- Annealed and Equilibrated Mix- 
ture protocols to generate the obstacle positions. As it 
has been shown in previous works, this theoretical scheme 
leads to an accurate description of the slow dynamics 
phenomena of fluids under different conditions. As can 
be seen in this work, SCGLE approach is capable to pre- 
dict the arrest transition diagrams in this kind of systems, 
independently of the protocol followed to generete obsta- 
cle positions. The required input are the static structure 
factors of the mixture which could be calculated using 
ROZ equations with an apropriate closure relation, once 
the protocol is defined. For EM systems, SCGLE arrest 
transition line is in excellent qualitative agreement with 
computer simulation results, showing a notorious reen- 
trant pocket which is well understood and has already 
been discussed in this work and in Ref. [Q] . On the other 
hand, for QA mixtures, we could compare our predictions 
with both computer simulations and another theoretical 
approach, Replica Mode Coupling Theory. The quali- 
tative agreement with computer simulation results is as 
good as in EM systems, but unlike RMCT predictions, 
the reentrant pocket is barely noticeable. In this way, 
we think that SCGLE has provided us with a theoretical 
approach to predict dynamic arrest in mobile-immobile 
mixtures independently of the structure of the obstacle 
matrix in a unified way, just providing the appropriate 
static inputs. 

This work was supported by Consejo Nacional de Cien- 
cia y Tecnologia through grants CB-20I0-C0I-156423 
and Red Nacional de Nanociencias y Nanotecnologia. 
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